Alkali Gases in Optical Lattices: A New Type of Quantum Crystals? 
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Similarities between alkali gases in optical lattices with non-integer occupation of the lattice 
sites and quantum crystals are explored. The analogy with the vacancy liquid (VL) provides an 
alternative explanation to the Mott transition for the recent experiment on the phase transition 
in the lattice. The VL can undergo BEC with T c within experimental reach. Direct and vacancy- 
assisted mechanisms of the band motion for hyperfine impurities are discussed. Large concentration 
of vacancies can result in the spatial decomposition of the system into pure hyperfine components. 
Below the vacancy condensation the impurity component resembles 3 He in 3 He — Hell mixtures. 



Recently, after the spectacular experimental discov- 
ery of Bose condensation, the study of alkali gases in 
traps has become the focal point in atomic, low temper- 
ature, and condensed matter physics. One of the most 
fascinating features is the possibility to see in an exper- 
iment some of phenomena that have been discussed ear- 
lier only within theoretical models (see review pj). An 
additional attraction is that the phenomena in ultracold 
alkali gases are incredibly rich and combine features in- 
herent to diverse condense matter and low temperature 
systems (Refs. and references therein). For example, 
BEC condensation in trapped gases resembles, but is not 
quite the same as the transition in other superfluid or 
superconducting systems [2j . Another example is the dy- 
namics of the hyperfine components which resembles the 
spin dynamics of spin-polarized quantum gases 

A new example is an ultracold alkali gas in an opti- 
cal lattice. Alkali atoms are almost localized in micro- 
scopically periodic potential wells induced by the Stark 
effect of interfering laser beams (Refs. and references 
therein) . The tunneling probability t between the wells is 
determined by the depth and size of the wells, i.e., by the 
intensity and the wavelength A of the beams. Since this 
intensity is adjustable, the atoms can be studied in a wide 
range of the tunneling frequencies t and effective masses 
in* ~ h 2 /ta 2 , from an almost free gas with a periodic 
perturbation to a well- localized "solid" (a = A/2 = ir/k 
is the lattice period). Another parameter is the on-site 
repulsion U for atoms inside the same well. 

The standard Hubbard model for electrons predicts Q 
that the alkali atoms in the optical lattice should ex- 
hibit an analog of the Mott metal-insulator transition at 
U/tPa 5.8z (z is the number of the nearest neighbors). At 
U 3> t, the system should become an "insulator" with- 
out interwell transitions (the particle tunneling between 
the lattices site increases the on-site energy by U and 
is energetically prohibitive). At U <C t, the on-site in- 
teraction does not restrict tunneling between the sites 
and the atoms are in the "metal" phase. Then, at suf- 
ficiently low temperature, the system can undergo BEC 
into a lattice superfluid. The energy parameters U, t, 
and the potential well Vq are often measured in units of 
the recoil energy E r = h 2 k 2 /2m. A typical example is 
@ t/E T ~ 0.07 and U/E r ~ 0.15 for V /E r = 15. 

This type of transition has been reported in Ref. |6J for 



lattices of the size a = n/k ~ 426 nm with recoil energy 
E r = h 2 k 2 /2m ~ 1 kHz. At low beam intensity, i.e., at 
large t, the experiment revealed the condensate peak in 
the center of the trap. This peak disappeared at small t 
(at Vo/E r between 13 and 22) which might indicate the 
transition to the Mott insulator (MI) phase. However, 
the identification of the high-C/ phase as the MI is not 
unambiguous. The MI can be observed only when the 
average number of atoms on the same lattice site is inte- 
ger. If the average population is fractional, the highest 
on-site energy states are not fully occupied. The tunnel- 
ing of the "excessive" particles from the site on which the 
highest level is occupied to the site with an unoccupied 
level cannot be banned by the on-site interaction. The 
tunneling of the "excessive" particle from site ri to the 
empty neighboring site r"2 increases the on-site energy by 
U on the site r2 while simultaneously decreasing it by U 
on the vacated site ri . Since both sites are translationally 
equivalent, this opens the way to the band motion of the 
" excessive" particles and to the existence of the partially 
filled conduction band. Then the lattice with non-integer 
occupation stays in the " metal" or " semiconductor" state 
even at large U with the " excessive" particles in the con- 
duction band. Ref. [|J contains the experimental proof 
of a large gap between the filled and conduction bands. 
However, it is difficult to conclude whether in equilibrium 
the conduction gap is empty or not. Below we suggest an 
alternative interpretation for Ref. based not on the 
analogy with the Mott transition, but on the analogy 
with the quantum crystals (QC). 

There is a strong similarity between the ultracold par- 
ticles in optical lattices and atoms in QC, such as solid 
helium, in which the tunneling is sufficiently high to en- 
sure band motion of atoms unless prohibited by large 
on-site repulsion (see review and references therein). 
In helium crystals, the atomic band motion is banned, 
as for all MI, when all the lattice sites are occupied by 
identical particles with the occupancy equal to one. If, 
however, some of the lattice sites are empty, nothing pro- 
hibits tunneling of atoms from the occupied onto the va- 
cant sites leading, as a result of translational symmetry, 
to the band motion of vacancies, i.e., to the formation 
of peculiar band quasiparticles - vacancy waves. Similar 
quasiparticles are formed when some of the atoms oc- 
cupy interstitial sites and can tunnel through the QC. 
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A slightly different situation occurs when some of the 
lattice sites are occupied by atoms of a different kind 
from the host matrix (impurities). The impurities can 
also tunnel through the QC despite the fact that each 
site still has the occupancy equal to one. The tunneling 
constant for impurities is smaller than for the vacancies 
since the exchange of places between the impurity and 
host atoms involves high-energy intermediate states with 
either double on-site occupancy or the atom in an inter- 
stitial position. The impurity-host exchanges could be so 
low that a more efficient mechanism of impurity motion 
could be the vacancy-assisted diffusion. The behavior of 
vacancy and impurity waves in QC is well understood 
0- However, the most exciting possibility in QC - su- 
perfluidity and BEC in the system of vacancy waves - 
has not been realized for "classical" QC, namely, solid 
4 He despite two decades of intensive efforts (Refs. 0, H 
and references therein). The reason is that there are no 
zero-temperature vacancies in solid He: with lowering 
temperature, the concentration of vacancies drops expo- 
nentially always remaining insufficient for BEC. 

The alkali atoms in optical lattices resemble QC with 
a very appealing difference: the BEC in the system of va- 
cancy or "impurity" waves could be within reach. When 
the occupancy of the individual wells is close to an in- 
teger, the system resembles QC with either a small con- 
centration of vacancies or "excessive" atoms. If the occu- 
pancy is slightly below the integer K, K = 1 + Int [n/N] , 
where n and iV are the densities of atoms and lattice sites, 
then the density of vacancies is n v = KN — n <C N. 
If the occupancy slightly exceeds the integer number, 
then the density of "excessive" atoms n e = N — n v is 
n e = n — Nlxit [n/N] <C N. Since the tunneling probabil- 
ities t are the same for vacancies and "excessive" atoms, 
t v = t e (in both cases, an atom tunnels to an empty 
site) , many properties of the system are symmetric with 
respect to the vacancies and excessive atoms. [This is not 
so for usual QC in which the lattice potential is built of 
the atom interaction and the potential relief is different 
for a vacancy and an "excessive" atom]. 

Below we consider the situation with large on-site in- 
teraction U when the lattice system with the integer site 
occupation would become a MI making the BEC tran- 
sition impossible. The analogy with QC automatically 
excludes fully occupied lowest on-site states and does not 
require a concept of counter-superfluidity [Tof. 

In the tight binding approximation for vacancies in a 
simple cubic lattice, 

e v (p) = A/2 - 2t v cos faa/h) (1) 

where A = 12t v is the bandwidth. At large U, the fixed 
chemical potential \x v is finite in contrast to fi v = for 
thermally-activated vacancies in solid helium. When n v 
is small, the BEC transition temperature for the vacan- 
cies can be determined using the standard equations for 
lattice gases with low band filling: 

T c = 6.6a 2 t v n 2 J 3 (2) 



and similar for excessive atoms. A good extrapolation 
between these two limiting cases is 

T c = 6.6a 4 t v n 2 J 3 (N - n v ) 2/3 . (3) 

With the above values of E r and a, the estimate for the 
BEC transition in the vacancy liquid (VL) is T c ~ 3 x 
I0~ 7 (t v /E r ) a;?/ 3 (1 — x v ) 2 ^ 3 K where x v = a 3 n v is the 
fraction of unoccupied highest on-site states. 

Though T c for the VL can be quite high, there are two 
reasons why T c is lower than the BEC temperature for 
a free gas. First, the density of participating particles is 
lower (only the vacancies or the excessive atoms in the 
highest on-site state are subject to the condensation) . In 
experiment [|| with the occupancy between 2 and 3, this 
leads, at least, to a factor 5~ 2 / 3 in T c and even stronger 
lowering of T c if the system is close to an integer occu- 
pancy. Second, the effective mass of vacancy waves or 
excessive particles, m* = h 2 /2ta 2 — m (E r /w 2 t) , could 
be much larger than the mass of the free atoms m. In 
experiment, one has limited control over the vacancy con- 
centration. On the other hand, the tunneling frequency 
depends exponentially on the intensity of the laser beams. 
This makes m* a readily adjustable parameter that can 
make the superfluid transition in the VL © observable. 

All this suggests a probable alternative to the Mott 
transition for experiment 6). At large t, the experiment 
confirmed the presence of BEC, probably in the "free" al- 
kali gas rather than in the VL. At small t, the experiment 
showed the absence of the condensate. However, the ex- 
periment, by design, cannot distinguish between the MI 
and the VL. It is very likely that the experiment showed 
the superfluid - VL transition rather than the superfluid 
- MI transition. 

The analogy with QC allows one make several other 
predictions. First, the role of "impurity waves" can be 
played by atoms in the different hyperfine states. Experi- 
mentally, such impurities can be studied by means similar 
to the NMR methods for 3 He diffusion in solid 4 He. In 
principle, impurities become band particles spread across 
the system with the tunneling frequency ti. However, 
when the on-site interaction U is large, i, for direct ex- 
change of impurities with the host atoms and, therefore, 
the impurity wave bandwidth become very small, of the 
order of ti ~ t 2 /U <C t v . When U is sufficiently large, 
such direct tunneling can be often disregarded. 

If the number of upper-state vacancies is noticeable, 
the vacancy-assisted processes dominate the impurity 
motion with an effective intersite tunneling rate ti ~ 
t v x v 3> t 2 /U (in this context, the asymmetry of the 
vacancy-assisted motion is not important). At T > 
T c , the vacancy-assisted processes are responsible for in- 
dependent tunneling transitions between the adjacent 
sites. This is not a band motion but a more traditional 
impurity diffusion with an effective diffusion coefficient 
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For the vacancy-assisted impurity tunneling, in contrast 
to the pure system (|3J|, and the mean free paths are 
not symmetric with respect to n v — > and n v — > N. 
In the former limit, the free paths are atomic while in 
the latter limit the impurities recover their band prop- 
erties with the large mean free path determined by the 
impurity-impurity scattering or the scattering by the few 
remaining upper-state host atoms. 

The situation changes dramatically after the vacancy 
system undergoes the superfluid transition © . Then the 
impurity becomes a completely delocalized quasiparticle 
in the vacancy superfluid background similar to 3 He im- 
purities in the superfluid 4 He [l2j. The effective mass 
of such quasiparticles at T — is m* ~ ft 2 /2t v x v and 
goes up with temperature with a decrease in density of 
the vacancy condensate. The interaction effects in this 
quasiparticle gas are negligible and the properties of the 
system can be evaluated using the standard equations 
for an ideal lattice gas of quasiparticles. At low enough 
temperatures, this impurity component of density rt{ will 
also undergo its own BEC with 



a 5 2/3 3 1/3 2/3™ 

b.o2t v an v n i ' = an v ' n i l c 



(5) 



where T c is the temperature for the vacancy condensa- 
tion (J2J. The emerging two-condensate system should 
exhibit properties similar to those of liquid 3 He — A He 
mixtures with two condensates below the 3 He transition 
[l3j| . Since this BEC in is based on the vacancy-assisted 
tunneling, this two-condensate system is different from 
the one considered in Ref. [To| . 

This picture of the vacancy-assisted impurity motion 
works well when the concentration of the hyperfine im- 
purities X{ = a 3 rii is low. At higher Xi the vacancy mo- 
tion in this translationally inhomogeneous environment 
is accompanied by the host-impurity permutations sup- 
pressing the band motion. This is similar to the vacancy 
motion in solid 3 He with disordered spin system. Then 
the vacancies autolocalize within the homogeneous do- 
mains of the size 



R = 



nl/5 



2m*NT [(x, - 1) In (1 - x t ) - x l Irixi] 



(6) 



which are filled by particles in one hyperfine state (Na- 
gaoka polarons). If the density of vacancies is large, 

1 /3 

n v R > 1, this should lead to the decomposition of the 
system into macroscopic hyperfine domains. The differ- 
ence between this decomposition and the vacancy-driven 
spin polarization of solid 3 He 0] is that the transition 
takes place when the concentration of the zero-point va- 
cancies and "polarization" (the concentration of hyper- 
fine components) are fixed. In contrast to the formation 
of dynamic, transient domains in experiment [3J, this de- 
composition leads to stationary domains. If the hyperfine 
impurities are bosons, this decomposition is not always 
necessary below the vacancy BEC. 

One feature of the optical lattices is quite different from 
more "traditional" QC such as helium. Since the periodic 



potential in QC is built of atomic interaction, the vacancy 
motion is tied to the deformation of the lattice. For the 
alkali atoms in optical lattices, the lattice is the external 
potential of the laser beams and the particle displacement 
is decoupled from the deformation of the lattice. As a re- 
sult, the low-frequency collective modes above and below 
BEC are decoupled from the lattice variables. 

An important issue for QC is the sensitivity of the 
narrow-band particles to external fields. Since the en- 
ergy of band particles cannot change by more than the 
bandwidth A, the external field f2 (r) makes the motion 
finite and localizes the particles in the area of the size 
Sr ~ aa, <r = A/ (dil/dr)a. For a usual QC, the im- 
portant fields are the particle interaction, lattice defor- 
mation, and external forces. In the optical lattice, the 
most important field is the trapping potential Q — \ctr 2 . 
In wide traps (small a), this trapping potential does 
not cause noticeable Umklapp processes and the overall 
Hamiltonian of a particle (or a vacancy) in the optical 
lattice, H — e (p) + |ar 2 , can be treated quasiclassi- 
cally. Analysis of this Hamiltonian can be performed in 
the momentum representation in which \ar 2 — ^h 2 aSI 2 y 
and the problem reduces to that for a particle with 
"mass" 1/a in the "potential" e (p). The quantum prob- 
lem is the simplest near the band minima where the 

quantized motion is harmonic with the characteristic fre- 
1/2 

quency u>q = (2tay / a/h. Since the effective (tunneling) 
mass to* = ti 2 /2ta 2 is larger than the mass of free par- 
ticles, this frequency is (to/to*) 1 ^ 2 times lower than for 
the free particles. The quantization in the trapping po- 
tential f2 (r) is usually not important and the motion is 
close to classical. At large a the motion is unrestricted. 
When a — ► 1, even the classical motion becomes com- 
pressed towards the multiwell shells around the center of 
the trap with a giving the thickness of the shell, to which 
the particle motion is restricted, in terms of the well size 
a. In experiment 0, ~ 75 (t/Er) 1 ^ 2 Hz and is small 
while er > 100 (t/E r ). The shells narrow to a single well 
layer at the beam intensity for which t/E r < 0.01. 

The inhomogeneity of the trap also leads to the non- 
uniform spatial redistribution of particles p|. If the 
change of the trapping potential from well to well is 
large in comparison with temperature (A/Tc 3> 1), the 
shells with the lower energy are fully filled, have integer 
population, and become the MI. The rest of the shells 
(most likely, but not necessarily, the outer ones) will have 
the non-integer population and resemble the VL with a 
rather large density of vacancies. In experiment [||, pa- 
rameter A/Ta ~ 10" 10 /T (with T in K) seems to be 
small meaning insignificant redistribution of particles be- 
tween the shells. Even if this parameter were large, the 
system would represent a thick shell with a > 100 (t/E r ) 
of coupled well layers in the quasi-2D VL state with the 
rest of the shells in the MI state with filled upper levels. 
This may actually increase the BEC temperature for the 
VL since the VL, though restricted to a lower number of 
shells, can have a higher density of vacancies. 

In summary, we explored the analogy between the al- 
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kali gases in optical lattices with non-integer occupation 
and large on-site interaction with QC. This analogy pro- 
vides an alternative explanation for experiment ja| as a 
transition between the BEC and VL states. The BEC 
transition for the VL is predicted. The transition temper- 
ature seems to be within experimental reach. The pres- 
ence of a large number of unoccupied states provides a 
vacancy-assisted mechanism for diffusion of hyperfine im- 
purities and can, sometimes, lead to a spatial decomposi- 



tion of the system into pure hyperfine components. The 
properties of the hyperfine mixture depend on whether 
the system is above or below BEC temperature for the 
VL. At even lower temperature one can observe the tran- 
sition to the state with two - vacancy and impurity - 
condensates which is different from Ref. [Tof . 

I am grateful to the referee for pointing out the im- 
portance of redistribution of particles in traps. The work 
was supported by NSF grant DMR-0077266. 
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